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Recently, a spatially extended excess of gamma rays collected by the Ferraz-LAT from the inner region of the 
Milky Way has been detected by different groups and with increasingly sophisticated techniques. Yet, any final 
conclusion about the morphology and spectral properties of such an extended diffuse emission are subject to a 
number of potentially critical uncertainties, related to the high density of cosmic rays, gas, magnetic fields and 
abundance of point sources. We will present a thorough study of the systematic uncertainties related to the 
modelling of diffuse background and to the propagation of cosmic rays in the inner part of our Galaxy. We will 
test a large set of models for the Galactic diffuse emission, generated by varying the propagation parameters 
within extreme conditions. By using those models in the fit of Fermz-LAT data as Galactic foreground, we 
will show that the gamma-ray excess survives and we will quantify the uncertainties on the excess emission 
morphology and energy spectrum. 


1. Introduction 


One of the most challenging results for indi¬ 
rect dark matter searches in recent years is the 
discovery of an excess emission in the gamma- 
ray flux from the center of our Galaxy. The 
first in dications o f suc h an ex cess d ate back _to 
2009 [Goodenough and Hooper! I2009L I Vitale et al.1 
2009}. Since then, several analyses of gamma- 
ray data from the Large A rea Telescop e aboard 
the Fermi satellite jGehrels and Michelsonl Il999l |. 
hereafter Fermi- LAT, claimed the existence of the 
excess above the sta ndard astrop hysical back ground 


at GeV energies Goodenough and Hooper 

2009, 

Hooper and Goodenough 

2011. Bovarskv et al. 

2011, 

Hooper and Linden 2011 

Abazaiian and Kaplinghat 

2012 

Macias and Gordon 2014, 

Abazaiian et al. 

2014 

, Davlan et al. 2014, Zhou et al. 

2013. The ex- 


cess emission results from anal yses o f bot h the inner 


few degrees of the Galaxv Abazaiian and Kaplinghat 

20121 Macias and Gordon! 

2014. Abazaiian et al. 

2014 Davlan et al. 2011 

, Zhou et al.1 2014. 

Gordon and Macias 20131 

and hipher lati- 

tudes Davlan et al. 2014 Hooper and Slatver 

20131 Huang et al. 2013], extending up to tens of 


degrees. Intriguingly, the observed spectral en¬ 
ergy distribution and the spatial properties of the 
Fermi GeV excess match the expectation for a signal 
from dark matter particles annihilating in the halo of 
the Milky Way. Nevertheless, some discussion about 
astrophysical explanations were put forward, as, for 
example, about the emission from a population of 
point-like sources below the telescope’s detection 
threshold [Hooper et al.1 120131 ICalore et al.l l2014al 
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ICholis et al.l 12014 iPetrovic et al.l 12014a |, or violent 
burst events at the Galactic center with injection 
of lep tons and/or protons s ome ki lo-/mega-years 
ago [Carlson and Profurnol I20l4 IPetrovic et al.l 
l2014bl| . 


Regardless of the possible interpretations, all anal¬ 
yses agree on the fact that an extra-emission over the 
standard astrophysical background is present in the 
inner region of the Galaxy. We stress here that the 
Galactic center is one of the most promising targets 
for dark matter searches since there the typically pre¬ 
dicted profiles for the dark matter distribution lead to 
the largest photon flux from dark matter origin. How¬ 
ever, the Galactic center is maybe the most challeng¬ 
ing target for dark matter searches: our knowledge of 
the conditions at the Galactic center is indeed very 
poor and the astrophysical background (from point 
sources as well as from diffuse emission processes) is 
affected by large uncertainties. 

A critical point is to answer the question “An ex¬ 
cess above what?”. The excess emission is defined 
with respect to specific astrophysical foregrounds and 
backgrounds, like the Galactic diffuse emission (which 
originates from the interactions of cosmic rays with 
gas and photons in the Galaxy), point-like and ex¬ 
tended sources. Those components should be mod¬ 
elled independently. Therefore, it is crucial to explore 
different foreground and background models in order 
to robustly identify and characterise the excess emis¬ 
sion. 


We will present here part of the analysis performed 
in ICalore et al.l 2014bj| . where we showed for the first 
time that the excess is statistically robust against 
theoretical model systematics, bracketed by exploring 
previously neglected uncertainties on the Galactic dif¬ 
fuse emission, and that the proper treatment of back¬ 
ground modelling uncertaint ies allows more free dom 
for models fitting the excess [Calore et al.ll2014cl ]. 
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Figure 1: Left panel: ICS and 7r°+Bremss contributions to the P6V11 background model. Right panel: Spectra 
predicted by a typical Galactic diffuse model for ICS, n° and bremsstrahlung emission. Fluxes in the 40° x 40° ROI, 
| 6 | > 2 °. 


2. On the importance of foreground 
modelling 

The dominant source of background for the Galac¬ 
tic center analysis is the emission originating from the 
interaction of cosmic rays with dust and gas in the 
Galaxy. The three main production mechanisms of 
Galactic diffuse gamma rays are: the Inverse Compton 
scattering (ICS) of electrons on low-energy ambient 
photons, the decay of abundantly produced neutral 
pions and the bremsstrahlung of electrons in the in¬ 
terstellar medium. Most of previous analyses adopted 
the same background model to describe the Galactic 
diffuse emission, namely the P6V11 background model, 
provided by the Fermi- LAT Collaboration for the sole 
purpose of point source analysis. 1 Using this model 
for analysis of extended sources introduces systematic 
effects that might lead to biased statements about the 
spectrum and morphology of the Fermi GeV excess 
emission. 

To visualise this effect, we decomposed the P6V11 
model in the main contributions to the Galactic 
diffuse emission. The spectra for ICS, n° and 
bremsstrahlung (that we consider as a unique com¬ 
ponent “7r°+Bremss”) are predicted by a standard 
mo del for cosm i c-ray p ropagation in the Galaxy 
Csee lCalore et al.1 |2014bj ] for more details). We fitted 
simultaneously the ICS and 7r°+Bremss components 
to P6V11 mock-data. From Figure [l] left panel, the 


J http://fermi.gsfc.nasa.gov/ssc/data/p6vll/access 
/lat/r ing_f or_FSSC_f inal4. pdf 


reader can see that an extremely hard ICS emission 
at energies >10 GeV is an intrinsic property of the 
P6V11 background model. The effect on any analysis 
that employs it as Galactic diffuse emission model is 
to over-subtract the ICS component at high energies, 
forcing the GeV excess spectrum to fall-off at > 10 
GeV. 

This exercise demonstrates the relevance of mod¬ 
elling separately the different contributions to the 
Galactic diffuse emission. Indeed, ICS, 7r° and 
bremsstrahlung possess intrinsically different mor¬ 
phologies because of the different targets that orig¬ 
inate these components: the gas for the 7r° and 
bremsstrahlung, and the interstellar radiation field 
for the ICS. Moreover, given the different cosmic-ray 
species responsible of the gamma-ray emission (pro¬ 
tons for 7T° and electrons for ICS and bremsstrahlung), 
also the way in which the morphology changes with 
energy is different for the three contributions. 

Such arguments strongly motivated the study of the 
variation of the spectral and morphological properties 
of the excess due to the modelling of the Galactic dif¬ 
fuse emission. 


3. Home-brew Galactic diffuse emission 

In order to robustly identify the excess despite of 
large variations in the foreground emission, we built 
a set of Galactic diffuse models by varying cosmic- 
ray propagation parameters within a given set of as¬ 
sumptions. We note that the observed emission re¬ 
sults from a line of sight integral and, as such, it re¬ 
ceives contributions from all distances. In particular, 
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Figure 2: Left panels: Count maps (the |6| > 2° cut and the point-source mask are clearly visible). Central panels: 
Residuals when subtracting all emission model templates. Right panels: Residuals when re-adding the GeV excess 
template absorbed emission. 


the emission that comes from the Galactic center is 
in the models we are adopting - relatively subdom¬ 
inant (about 10%) and the Galactic diffuse emission 
is dominated by local processes. Therefore, our work 
should read as the characterisation of the uncertainties 
due to the Galactic gamma-ray emissivity along the 
line of sight. We worked in the framework of steady 
state solutions to the transport equation of cosmic- 
ray propagation in the Galaxy. Homogeneous diffu¬ 
sion, re-acceleration and convection were con si dered. 
We ad opted models from the set of lAckermann et al.1 
[2012j] to test variations of the diffusion zone geome¬ 
try, the source distribution, the spin temperature and 
the magnitude cut (for an explanation of cosmic-ray 
pro pagation parameters and their range of variation 
see iGalore et al.l |2014bj |). Additionally, we generated 
our own Galactic diffuse models using Galprop v54 
(webrun version). With those models, we explored 
the remaining uncertainties related to the diffusion 
coefficient, re-acceleration, convection, interstellar ra¬ 
diation field, and Galactic center magnetic field dis¬ 
tributions. In total, we built a set of about 60 models 
for the Galactic diffuse emission that test “extreme” 
variations in the parameter space. We here quote the 
explored parameter ranges: 

• geometry of the diffusion zone: 4 < zd < 10 kpc 
and rp = 20 or 30 kpc; 

• source distributions: SNR, pulsars, OB stars; 

• diffusion coefficient at 4 GV: Dq = 2 — 60 x 10 28 

2 —1 
cm s ; 


• interstellar radiation field model factors (for op¬ 
tical and infrared emission): 0.5 - 1.5; 

• magnetic field parameters: 5 < r c < 10 kpc, 1 < 
z c < 2 kpc, and 5.8 < B(r = 0,0 = 0) < 117 
/rG. 

We note that we did not test those models against 
local cosmic-ray data and large scale diffuse gamma- 
ray data (or even microwave data). 

As already mentioned, we made a few simplify¬ 
ing assumptions that we summarise below and that 
will become relevant for future refined analyses of the 
GeV excess : (i) homogeneity and isotropy of cosmic- 
ray diffusion, re-acceleration, and convection; (ii) ra¬ 
dial symmetry of cosmic-ray source distribution in the 
Galactic disk (i.e. no modelling of the spiral arms), 
and same source distribution for different cosmic-ray- 
species sources; (iii) steady state regime, excluding 
transient phenomena as, for example, burst events. 


4. The data analysis 


In order to analyse gamma rays collected by the 
Fermi -LAT from the inner Galaxy, we adopt a 
template-based multi-linear r egress io n technique, see , 


template-based multi-linear regression tccnmquc , see , 
for example, iDobler et al.l 201flj . ISu et all [2010j |. 
The data sample corresponds to 284 weeks of repro¬ 
cessed Fermi -LAT data (from 4 August 2008 on) in 
the energy range 300 MeV - 500 GeV. The Region- 
Of-Interest (ROI), i.e. the inner Galaxy, is defined as 


\£\ < 20° and 2° < |6| < 20° 


(1) 


• Alfven speed: va = 0 — 100 kms 1 ; 

• gradient of convection velocity: dv/dz = 0 - 500 
kins -1 kpc -1 ; 


The choice of the latitude cut is such to avoid the large 
contamination of point sources in the innermost few 
degrees, where the source confusion is very high. We 
prepare the data according to standard prescriptions 
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provided by the Fermi-Science-Support-Center. We 
binned the data on an healpix grid with resolution 
parameter n s i ze = 256, for each energy bin (24 in to¬ 
tal) defined in a way such to guarantee good statistics 
also at the highest energies. 

We compare the data maps with the model maps, 
obtained by the superposition of the different tem¬ 
plates adopted in the analysis (see below). The 
best-fit normalisation of each model template is de¬ 
rived through a maximum likelihood method, based 
on t he P oisson likelihood function (cf. Eq. (2.3) 
in ICalore et al.1 [2014bj |h 

The spatial model templates adopted in the analysis 
are: 


Point-like so urces template as derived from the 
2FGL A in lo et ahl 2011 1, with fixed spectra and 
flux normalisations. 


Fermi bubbles modelled by a uniform¬ 
brightness spati al template with bubbles’ 


edges as in ISu et al.l 201 0j. 


• Isotropic gamma-ray diffuse background with 
uniform-brightness emission template. 


• Galactic diffuse emission ICS and 7 To+Bremss 
independent templates as modelled from Sec. EH 


• GeV excess template whose volume emissivity is 
parametrised by the spherically symmetric gen¬ 
eralized NFW profile, 


P(r) = Ps 


(r/r s ) 7 
(1 -(- r/r s ) 3 ~'i ’ 


( 2 ) 


squared, and with (best-fit) spectral index 7 = 
1.2. This choice is clearly motivated by the dark 
matter annihilation interpretation of the GeV 
excess, although we tested a large range of vari¬ 
ation for the profile parameters. 


The fitted spectra of the Fermi bubbles and isotropic 
diffuse background templates are constrained to vary 
wit hin the measured s pectra from iFranckowiakl (2013] 
and I Ackermannl 2012|] . respectively. 

In the analysis, we introduced the following techni¬ 
cal improvements: a non-logarithmic energy binning 
such to counterbalance the reduced photon statistics 
above 10 GeV, a weighted adaptive masking of point 
sources, and the full treatment of the Fermi -LAT 
point spread function. 


5. Selection of main results 

In this section, we present a selection of t he res ults 
of the a nalysis, and we refer the reader to lCalore et al.l 
2014b| ] for a thorough explanation of our findings. 
Figure [2] represents the residual (i.e. data - model 
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Figure 3: Spectral energy distribution of the GeV excess 
template, for a generalised NFW profile with an inner 
slope 7 = 1.2. The yellow band corresponds to all of the 
60 GDE models. Two models are highlighted: the model 
that provides the best fit to the data (model F, green 
points ) and a reference model (red points ). 


counts) emission obtained when subtracting from the 
raw data the emission associated with the model tem¬ 
plates (central panel). The residuals are at the level 
of 20% in the whole ROI, but, when the GeV excess 
template associated to the model is re-added (right 
panel), the residuals in the central region of the ROI 
increase significantly, attesting the presence of the ex¬ 
cess, which is, after the other components are sub¬ 
tracted, the most pronounced large-scale excess in our 
ROI. 

Figure [5] represents the spectral energy distribution 
of the excess emission, i.e. the emission absorbed by 
the GeV excess template during the fitting procedure. 
The yellow band results from all the adopted Galactic 
diffuse models. Such a band brackets the uncertainty 
due to the theoretical modelling of the Galactic dif¬ 
fuse emission and affecting the extraction of the GeV 
excess spectrum. The GeV excess emission is found 
to be remarkably stable against the tested variations 
of the Galactic foreground. The typical GeV excess 
spectrum shows a rising below 1 GeV (with a spec¬ 
tral index harder than ~ 2 for all Galactic diffuse 
models) and features a peak at energies around 1-3 
GeV. Despite previous analyses, at higher energies, 
the spectrum is described by a power-law with slope 
~ -2.6. The coloured data points indicate the spec¬ 
trum (with statistical errors) that corresponds to the 
best-fit Galactic diffuse model ( model F) and another 
exemplary model discussed in ICalore et al.l [2014bl ] 
(model A). 

The envelope of the yellow lines corresponds to the 
theoretical model uncertainty , which is due to the vari¬ 
ation induced by the Galactic diffuse modelling. Such 
uncertainty is, at all energies, larger than the statis¬ 
tical errors, indicating the importance of the proper 
treatment of background model systematics. 
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Figure 4: GeV excess emission spectrum, together with statistical and systematical errors, for model F (i.e. best-fit 
model). Several spectral models have been fitted to the data. All the spectra (except the t + t~) provide a quite 
reasonable fit to the data. This is due to the correlation of the systematic errors (see text). 


Spectrum 

Parameters 


X 2 /dof p- value 

broken PL 

ai = 1.42±°; 22 , a 2 = 2.63±°;J 3 5 , F break = 

2.06t°; 23 GeV 

1.06 

0.39 

DM XX bb 

(av) = 1.76±°' 2 ? x 1CT 26 cm 3 s _1 , m x = 

= 49+®'4 GeV 

1.08 

0.36 

DM XX cc 

(ov) = 1.25to 2 8 x 10 ” 26 cm 3 s -1 , m x = 

38.2±+® GeV 

1.07 

0.37 

PL with exp. cutoff 

E cut = 2.53tJ; 3 7 GeV, a = 0.945j)°; 36 

1.37 

0.12 

DM xx —1 t + t~ 

(av) = 0.337jlo'o48 x 10" 26 cm 3 s" 1 , m x = 

= 9.961^ GeV 

1.52 

0.06 


Table I Spectral fits to the GeV excess spectrum, with ±1 <t errors. We show best-ft parameters, reduced x 2 , and 
corresponding p-value. 


As it can be already deduced from the residual 
plots, the Galactic diffuse models tested in the present 
analysis do not describe the data at the statistical 
level, but, still, they show large residuals in the ROI. 
Indeed, although the reduced y 2 for the best-fit Galac¬ 
tic diffuse model (model F) in the energy range from 
500 MeV to 3.31 GeV is close to one (~ 1.10) be¬ 
cause of the large number of free parameters, the cor¬ 
responding p-value is ridiculously small, ~ 10~ 300 . 


On the base of this argument, it is important to find 
an alternative way of assessing t he systematics uncer - 
tainties affecting the excess. In ICalore et al.l 2014bj| . 
we relied on an empirical method to derive model sys¬ 
tematics due to how well the different Galactic diffuse 
models describe the data along the disk, away from 
the Galactic center. The derivation and definition of 
the empirical model systematics were presented dur¬ 
ing this conference in a complementary talk, “Robust 
Identification of the GeV Galactic Center Excess at 


Higher Latitudes”. 2 Quantifying the background em¬ 
pirical model systematics turned out to be crucial for 
making statistics based claims on the possible inter¬ 
pretations of the excess. 


6. Interpretations 

As explained in Sec. [I] several interpretations have 
been proposed, ranging from purely astrophysical to 
dark matter explanations. As a first constraint, the 
predicted model spectrum must provide a good fit to 
the GeV excess spectrum. We performed parametric 
fits to the GeV excess observed spectrum fully taking 
into account the systematic uncertainties. 


2 C. Weniger et al., proceedings RICAP-14 (to appear soon). 
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To do so, we made use of a \ 2 function with a non¬ 
diagonal covariance matrix, which models the corre¬ 
lated empirical model systematics. The y 2 function 
used writes as: 



with £- 1 the covariance matrix. The covariance con- 

l J 

tains model uncertainties that were derived from the 
size of typical residuals along the Galactic disk. They 
amount to variations in the excess template that are 
similar to the ones shown in Figure 0 and are illus¬ 
trated in Figure |4] by the yellow boxes. For details we 
refer the reader to iGalore et al.l 2014b ]. 

We tested several spectra that are related to the 
GeV excess viable interpretations. TablcUsummarises 
our findings. In particular, parametric fits with corre¬ 
lated errors show equal preference for a broken power- 
law spectrum and for the spectrum from dark matter 
annihilation into b-quarks. Remarkably, the p-value 
for a spectrum due to dark matter annihilation into 
r-leptons is higher the 0.05. The reason for which 
dark matter annihilation spectra provide good fit to 
the GeV excess is due to the fact that systematics er¬ 
rors are correlated in energy and can be understood in 
terms of t he covariance matrix (we refer the interested 
reader to ICalore et al.l |2014cj] h 


7. Conclusion 

The analysis of the Fermi- LAT data performed 
in ICalore et al.l [2014bj] confirmed the presence of an 
excess emission in the inner Galaxy and some of its, 
previously found, specific properties, such as the 2- 
3 GeV peaked spectral energy distribution, the ex¬ 
tension to high latitudes and the compatibility with 
a spherically symmetric spatial distribution. Those 
properties were demonstrated to be remarkably stable 
against theoretical model systematics, due to the vari¬ 
ations in the modelling of the Galactic diffuse emis¬ 
sion. 

We assess empirical model systematics from a scan 
of the gamma-ray flux along the disk and we used 
those uncertainties as a proxy for the systematics af¬ 
fecting the GeV excess at the Galactic center. 

Contrary to previous results, we do not confirm the 
fall-off of the GeV excess spectrum at E > 10 GeV, 
but we do find a high energy tail of the spectrum ex¬ 
tending up to 100 GeV. However, when we properly 
treat model systematics and include them in the spec¬ 
tral fits as correlated errors, we demonstrated that it 
is possible to equally well fit the excess spectrum with 
both a broken power-law and a gamma-ray spectrum 
typically expected from dark matter particles annihi¬ 
lation into bb final states. This implies a large, pre¬ 


viously neglected, freedom for models fitting the GeV 
excess. 
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